EPJ Web of Conferences will be set by the publisher 
DOI: will be set by the publisher 

© Owned by the authors, published by EDP Sciences, 2012 



Challenges in the modelling of tidal disruption events lightcurves 
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Abstract. In this contribution, I review the recent developments on the modelling of the lightcurve of tidal 
disruption events. Our understanding has evolved significantly from the earlier seminal results that imply a 
simple power-law decay of the bolometric light curve as t' 5 ' 3 . We now know that the details of the rise to 
the peak of the lightcurve is determined mainly by the internal structure of the disrupted star. We also have 
improved models for the disc thermal emission, showing that in this case the decline of the luminosity with 
time should be much flatter than the standard r 5/3 law, especially in optical and UV wavelengths, while the 
X-ray lightcurve is generally best suited to track the bolometric one. Finally, we are just starting to explore 
the interesting general relativistic effects that might arise for such events, for which the tidal radius lies very 
close to the black hole event horizon. In particular, I describe here some possible evidences for relativistic 
Lense-Thirring precession from the light curve of the event Swift J 1644. 



1 Introduction 

The theory of tidal disruption events (TDE) by supermas- 
sive black holes has been developed since the early '80s 
with the seminal work of Carter and Luminet HI 0] Lacy, 
Townes and Hollenbach [3], Rees Jj] and Phinney 
However, for several years little progress has been made, 
except for some important contribution from the numeri- 
cal point of view JaUl] and on the analysis of the accretion 
regime ]8[]. 

In recent years, flares from the nuclei of quiescent 
gala xies, detected either in X-rays |@- 11 ] and in the UV 
1 12l4l4ll have been interpreted as arising from TDEs, es- 
sentially by a simple comparison of the light curve of the 
event to the expected r 5 ^ law [4]. This has prompted sev- 
eral new theoretical investigations that aimed at clarifying 
the conditions under which one should really expect such 
a power-law decline in the luminosity of the event. 

The field has further expanded once it has been re- 
alised observationally that such events might produce 
powerful relativistic jets Il5lll6ll . whose modelling clearly 
needs to go beyond the simple estimate of the dynamics 
of the infalling debris and the disc thermal emission Inll . 
Additionally, the larger scale outflows that might be con- 
nected with the early super-Eddington fallback phase is 
expected to produce a characteristic spectroscopic signa- 
ture dO. 

In this framework, modelling the dynamics of the de- 
bris is still a key issue in order to understand the time evo- 
lution of the accretion rate onto the black hole, which is 
an essential input parameter for any detailed model of the 
system emission properties. In this contribution I will first 
describe the basic results and the standard theory behind it 



and I will then focus on a number of recent developments 
that have somewhat changed the basic expectations of the 
simplest models. 

2 Light curve modelling 
2.1 Classic results 

The tidal radius beyond which stellar self-gravity is not 
able to counteract the black hole's tidal field and keep the 
star together is defined as: 



1/3 



R, 



(1) 



where Mh is the black hole mass, M, the stellar mass and 
R t the stellar radius. Note that for typical parameters this is 
not much larger than the black hole's Schwarzschild radius 

^ * 23M" 2/3 , 
Rs 6 



(2) 
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where M 6 = M h /1O 6 M , which implies that black holes 
with masses above a 1O 8 M are not expected to produce 
any TDE. It should be noted however that the above limit 
is only valid in Newtonian dynamics. General relativistic 
calculations have shown 112 111 that highly spinning black 
holes are able to produce stellar disruption of solar type 
stars for masses well above 1O 8 M . This is a first example 
that shows how important general relativistic effects can 
be in describing TDEs, and thus be a powerful probe of 
black holes properties otherwise difficult to measure, such 
as their spin. 

Often, the orbit of the incoming star is described as a 
simple parabolic orbit, with a pericenter distance R p close 
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Figure 1. Left: fallback rate as a function of time for the disruption of three polytropic stars, with progressively decreasing values 
of 7 = 5/3 (solid), y = 1.4 (long-dashed), and y = 4/3 (short-dashed). Right: corresponding instantaneous value of the logarithmic 
derivative of M(t). More compressible, and centrally concentrated stars, have a more gentle rise to the peak, and the t' 5 " regime is 
reached later. From l2(5ll . 



to the tidal radius. One can define a "penetration factor" 
/3 = RJRp, such that when (3 < 1 tidal disruption should 
not occur. Since the details of the light curve essentially 
depend on the conditions at the tidal radius (see below), in 
the following we will only restrict to the case where /? = 1 
so that the tidal radius coincides with the pericenter of the 
stellar orbit. 

After disruption, the stellar debris are launched into 
very eccentric orbits (with 1 - 4R*/R p < e < 1) with 
a spread of orbital energies, and gradually return to peri- 
center, where they circularise and form an accretion disc 
at /? c j rc « 2R p . During this circularisation phase the de- 
bris lose a fair fraction of their energy, which might be 
in principle detectable. Note, however, that the ratio be- 
tween the luminosity released during circularisation and 
that released by the subsequent accretion of the debris is 
expected to be 



potential energy spread across a stellar radius: 



^ISCO 
Rc\xc 



0.06/3M1 13 , 



(3) 



and is thus generally neglected in most calculations. For 
large black hole masses, or for deeply penetrating events, 
this contribution might become significant. 

The fallback rate of the debris is easily derived from 
the the distribution of their orbital energies dM/dE and 
Kepler's third law (the effects of general relativity can also 
been included, I22I0 : 



(27rGM h ) 2/3 dM^ 
3 dE~' 



5/3 



(4) 



which, in the limit that the energy distribution of the de- 
bris is flat, gives the standard r 5 ^ 3 decay law, almost uni- 
versally associated with this type of events [4]. A simple 
order of magnitude estimate for dM/dE is obtained by as- 
suming that the debris are distributed uniformly within the 



dM 
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2.2 Recent developments 

2.2. 1 Effects of stellar structure 

The theory described above highlights the fundamental 
role played by the distribution of orbital energies of the 
debris. Clearly the assumption that it is flat is a first 
crude approximation, and the actual shape of the distribu- 
tion, and hence the actual time dependence of the fallback 
rate needs to be computed. The issue has been initially 
treated from a purely numerical point of view by Evans 
and Kochanek [6], who, somewhat surprisingly, appeared 
to confirm that indeed the hypothesis of a flat distribution 
originally made for simplicity by Rees Wl was correct. 
More recently, Lodato, King and Pringle B20I1 have revis- 
ited the problem both from an analytical and a numerical 
point of view. 

The analytical model is based on the assumption that 
the structure of the disrupted star is essentially unperturbed 
until it reaches the tidal radius. The assumption is clearly 
approximate, but is needed in order to compute analyti- 
cally the resulting energy distribution, which then needs 
to be compared to numerical simulations. The big advan- 
tage of such analytical model is that it essentially allows 
to compute the debris energy distribution, and hence the 
fallback rate, starting from the internal density profile of 
the unperturbed disrupted star: 

r>)/d/, (6) 

dE GM h X HK 

where r = ER^/GM^. The main result of this theory is that 
more compressible stars, which are more centrally concen- 
trated, have a more gentle rise to the peak, in the sense that 
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Figure 2. Expected fallback rates for the disruption of polytropic 
stars, with progressively increasing values of y = 1.4 (solid), 
y = 1.5 (short-dashed), y = 5/3 (dot-dashed) andy = 1.8 (long- 
dashed), scaled to solar values. More centrally concentrated stars 
(with a lower y) have a more gentle rise to the peak and peak at an 
earlier time with respect to more incompressible models. From 



120]. 



the logarithmic gradient of the light curve reaches the ex- 
pected -5/3 value at later times (see Fig. Q}. 

Lodato, King & Pringle lEoll have also tested this 
model numerically, by running SPH simulations of the dis- 
ruption of simple polytropic stars, with /3 — 1 . These nu- 
merical tests reproduce to a large extent the predictions of 
the model, hence confirming that the analytical model is 
a good first approximation. However, there are two im- 
portant changes to be taken into account: (a) the effective 
reduction of the stellar self-gravity as it approaches the 
tidal radius determines an inflation of the star itself and 
(b) shocks, that develop preferentially in more compress- 
ible stars, tend to push matter at the edges of the energy 
distribution, thus further modifying the fallback rate. The 
net result of these two effects is to delay the time of the 
peak for the less centrally concentrated stars (see Fig. 
Such results have also been confirmed by later investiga- 
tions of the same problem [23]. 

Form an observational point of view, the models above 
have been used to describe the UV emission of several 
TDE candidates lfl3ll . In particular, very recently Gezari et 
al have used the arguments above to conclude that a 
peculiar TDE flare can be attributed to the disruption of the 
helium core of a red giant, whose envelope had been previ- 
ously removed (possibly by the same tidal mechanism that 
later caused its complete disruption). All such analyses 
assume that the UV an optical light curve is proportional 
to the fallback rate. Whether such assumption is correct 
should be discussed in greater detail (see next section). 

It is worth noting that the above investigations have 
only considered the case where the pericenter of the orbit 
coincides with the tidal radius, and hence /3 — 1. How 
large is the energy spread of the debris for more penetrat- 



Figure 3. Bolometric, X-ray (0.2 keV), FUV and optical light 
curves of the tidal disruption of a solar mass star by a 10 6 M G 
black hole. Only the disc thermal emission is shown. The green 
lines indicate a r s/3 and a r s/n decline, respectively. From [25]. 



ing encounters? Is it the relatively small spread as deter- 
mined at the tidal radius or does it rather reflect the con- 
ditions at pericenter? The problem has been first investi- 
gated by Sari et al [24] in the related context of tidal dis- 
ruption of binary systems by a supermassive black hole, 
and later confirmed also in the context of stellar disrup- 
tion [23]. The conclusion is that the energy spread is de- 
termined by the conditions at the tidal radius, and is thus 
much smaller than what would be naively predicted if one 
computes it based on the conditions at the pericenter. 

2.2.2 Modelling the thermal disc emission 

Up to now, I have just considered the modelling of the 
fallback rate, which shows, at least after the peak of the 
emission, the typical r 5 ^ decline with time. How does 
this convert into luminosity? 

First of all, note that the fallback rate is the rate at 
which the debris return to pericenter and circularise, form- 
ing an accretion disc. As noted above, most of the lumi- 
nosity associated with a tidal disruption event is released 
as the debris accrete from the circularisation radius onto 
the black hole through an accretion disc. It is thus the 
accretion rate M onto the black hole that determines the 
luminosity of the flare. This in general is not equal to the 
fallback rate Mfb, unless the viscous timescale in the disc 
is much smaller than the timescale related to the fallback 
of the debris. We will discuss this condition more in detail 
in the following section. For the moment, let us just make 
the usual and simplifying assumption that M = Mfb- 

In this case, the fallback rate essentially sets the scale 
for the bolometric luminosity of the system, trhough 
L = T]Mf\,c 2 oc r 5/i , where r\ is the accretion efficiency. 
However, the light curves at specific wavelengths might 
strongly differ for the above behaviour, given that the sys- 
tem also cools down significantly while fading. The issue 
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Figure 4. The expected timescale for the transition from a 
fallback dominated accretion regime into a viscosity dominated 
regime, as a function of the black hole mass for various choices 
of a. 



has been studied in detail by Lodato & Rossi [25], who 
describe the time evolution of the thermal disc emission at 
various wavelengths. It can be thus demonstrated that the 
light curve, at any given wavelength, should go through 
three different phases. (1) Initially, when the disc is very 
hot, the observed wavelength might fall on the Rayleigh- 
Jeans tail of the blackbody spectrum produced by the disc. 
The luminosity L v in this case is thus proportional to the 
temperature L v oc Ta sc M^ 4 oc r 5 ^ 12 (note that this spe- 
cific result has also been previously suggested by IU8I0 . 
(2) As the disc cools, the observed wavelength might cor- 
respond roughly to the peak of the blackbody spectrum. It 
is in this phase that one should expect the single wave- 
length light curve to decline as r 5 ^, as the bolometric 
light curve. (3) Finally, once the observed wavelength 
moves to the Wien part of the spectrum, the light curve will 
start to rapidly fade off exponentially. The first of these 
three phases might be absent if the observed wavelength is 
high, for example in X-rays. These three behaviours can 
be clearly seen in Fig. [3] 

The above results thus indicate that X-ray observations 
are most likely to display the canonical f~ 5 ^ 3 decline, while 
the UV and optical light curve should display a much shal- 
lower light curve. Such a shallower decay in optical and 
UV has been observed in the case of Swift J2058 [16]. 



2.2.3 Fallback versus viscosity dominated accretion 

As mentioned above, one key assumption of the tradi- 
tional modelling of TDE is that the viscous timescale in 
the disc is much smaller than the return timescale of the 
debris. In this way, the disc can be described by a se- 
ries of steady state configurations with varying M, at each 
time equal to the incoming fallback rate Mr,. At the oppo- 
site end, some models |18[] have considered the behaviour 
of a system where the fallback material forms a narrow 
ring, which at later times starts accreting onto the black 
hole. Such models are appropriate only when the viscous 
timescale in the disc is always much longer than the fall- 
back timescale. 




Figure 5. Evolution of the instantaneous power law index of 
the light curve as a function of time for several choices of f tnffls 
(decreasing from top to bottom). The x-axis shows log(?/f m i„), 
where f min is the return time of the most tightly bound debris, 
marking the beginning of the flare. For the case indicated with a 
black line (explored here for completeness), f tnms is comparable 
to f min so that the t 5/3 plateau does not occur. 



In reality, neither model is correct. At earlier times, 
when the accretion rate is close to Eddington and the disc 
is hot and thick, it is expected that the viscous timescale is 
indeed much shorter than the fallback time, while at later 
times, once the accretion rate has fallen off significantly 
and the disc has cooled down, the ordering of the relevant 
timescales should be reversed, and the TDE should tran- 
sition from a. fallback dominated accretion regime, where 
the accretion rate coincides with Mq» into a viscosity dom- 
inated accretion regime, where it is the internal disc vis- 
cosity that determines the time evolution of M. 

The timescale for such transition can be easily com- 
puted based on the following argument. The timescale 
over which the fallback rate evolves is: 



Mo, 



5/3 



(7) 



On the other hand, the viscous timescale in the disc is 
proportional to (H/R)~ 2 (where H is the disc thickness). 
At the high accretion rates predicted for TDE the disc is 
likely to be radiation pressure dominated and thus H/R oc 
M/MEdA- We thus have: 



a\R) 



(8) 



where a is the disc viscosity coefficient and thus grows 
with time much faster than ffb. Thus, even if at small times 
t v <<c ffb, at later times we should expect a transition. The 
time at which the two timescales become equal is: 

/ a \ 3 / 5 / M h Y 1/5 ( MA 3 ' 2 ( RA~ 3 ' 2 

(9) 

which shows that for most cases this transition actually 
occurs after several decades from the disruption, once the 
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accretion rate has dropped too much to be significant (cf. 
Hill ). Figure 3] shows this timescale as a function of the 
black hole mass for several values of a. 

Detailed time-dependent calculations of this effect 
(Lodato, in preparation) show that significant departures 
(to a 20% level) to the canonical r 5 ^ decline can occur al- 
ready at a fraction of ? trans as small as 0.1. Figure [5] shows 
the results of one such calculations for various values of 
^trails- In particular we plot the instantaneous power-law 
index of the light curve as a function of time. One can 
see that after a plateau at -5/3, the light curve steepens 
significantly, once accretion through the disc slows down. 

3 Modelling Swift J1644: lightcurve and 
precession 

Swift J 1644 has been one of the most spectacular TDE 
reported to date, showing for the first time that powerful 



Swift/XRT dala o1 swill j1 64449.3+573 
blue: WT- red: PC 



at powi 

jets can be associated with such events H151 126H28I1 . A 



detailed modelling of the spectral energy distribution can 
reveal what are the main emission mechanisms for the jet 
emission lfl5ll . An analysis of the time evolution of the 
light curve can, on the other hand, give us clues on the dy- 
namical features of the event. Such a dynamical analysis 
has been attempted [29] based on the early light curve. 

We can attempt here a new analysis based on the light 
curve during the first half a year of activity. Figure[6]shows 
the XRT count rate (red and blue lines) for the event during 
the first six months. The system is very variable, but a 
simple t " can roughly reproduce the data (black line). 
This has been obtained by assuming that the jet luminosity 
scales with the fallback rate: 



Lx(f) = fM^ 2 



(10) 



where / is a scaling factor that might not be too differ- 
ent than unity, since while it is true that only a fraction of 
the accretion energy is expected to be released from the 
jet, one should also consider that the X-ray luminosity is 
boosted by relativistic beaming. We further assume the 
following behaviour of Mfb with time: 



Mfb(f) = M peak 1 + 



-5/3 



(ID 



where f m j n is the return time of the most bound debris and 
one should remember that in this case t — marks the 
beginning of the outburst (the BAT trigger) rather than the 
pericenter passage of the star, as before. 

This comparison allows us to estimate the following 
parameters: to, ~ 19 days, fM pe! ±c 2 ss 4 10 46 ergs/sec, 
M h * 2/ 10 6 M G , AT, ~ 0.1M o //. We can thus see that 
the parameters estimated in this way are in line with the 
expectations from such event, for an / factor not too dif- 
ferent from unity. 

Swift J 1644 might also show some additional interest- 
ing features connected with relativistic effects. The X-ray 
lightcurve indeed shows an interesting 2.7 days periodicity 
1 27] that might be connected with jet precession. Indeed, 
Lei et al 13011 have suggested that such modulation is linked 
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Figure 6. The blue and red curves show the XRT light curve 
of Swift J 1644, while the black line is a rough fit to it, using a 
simple f -5 ' 3 law. 




Figure 7. The spin parameter a of the black hole in Swift J 1644, 
as implied by assuming that the 2.7 days variability of the X-ray 
lightcurve is due to disc precession caused by the Lense-Thirring 
effect. Here p is the power-law index of the surface density pro- 
file (small p implying flat profiles). The values of the assumed 
black hole mass are indicated. 



to precession induced by the Bardeen-Petterson effect due 
to a warping of the disc caused by the torque associated 
with the spin of the accreting black hole. However, this is 
probably not the case because: (a) the Bardeen-Petterson 
effect does not lead to a steady precession. The inner disc 
rapidly gets aligned with the black hole spin, and the fi- 
nal realignment of the latter with the outer disc occurs 
on a timescale comparable with the precession timescale 
1 31]; (b) the Bardeen-Petterson effect occurs for very thin 
and viscous discs, which is not the case for the almost- 
Eddington accretion rates involved in TDE. 
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In reality, though, a steady precession of the disc is ex- 
pected in the case of thick configurations (see also 
where the warp propagates in the disc as a wave, provided 
that the disc radial extent is short enough (which is also 
expected in this case). The precession frequency is ob- 
tained simply by balancing the overall torque produced in 
the disc by the frame-dragging effect, and the disc inertia: 



J Q. LT (R)L(R)RdR 
J L(R)RdR 



(12) 



where Qlt is the Lense-Thirring precession rate [31]. 
Since Qlt depends on the black hole spin parameter a, one 
can use the known precession rate £2 p to infer the black 
hole spin. In particular, we can assume that the disc sur- 
fece density has a simple power-law dependence on radius 
Z,(R) oz R~P, where p is a free parameter. Figure [7] shows 
the implied value of a as a function of p for several values 
of the black hole mass. We thus see that, unless the black 
hole mass is small and the surface density profile is very 
flat, we estimate a rather low value for a ~ 0.2 - 0.4 with 
a relatively weak dependence on the free parameters. 

4 Conclusions 



The modelling of the light curve of tidal disruption events 
has moved significantly from the earlier standard results 
implying a r 5 ^ decline. We now know much better how 
to model the initial rise of the lightcurve, and how it is con- 
nected with the internal structure of the disrupted star. We 
begin to model in greater details the specific lightcurves at 
various wavelengths, although further work is required in 
this context to understand why in some cases the observed 
UV lightcurves appear not to cool down with time as ex- 
pected from a thermal disc emission. Finally we are just 
now starting to understand how relativistic effects (for ex- 
ample connected with the black hole spin) can affect the 
light curve of the system. In this respect, TDE offer a 
unique way to probe relativistic effects, given that the tidal 
radius is often very close to the black hole event horizon. 
Further specific modelling of such effects is going to be 
essential in order to correctly interpret upcoming observa- 
tions. 
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